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Based on the Pauli spin operators we develop the notion of the spin-correlation 
matrix for the two-qubit system. If this matrix is non-zero, the measure of the 
correlation between the qubits is the average of the non-zero elements. Trivially, for 
zero matrix the bipartite is uncorrelated. This criterion turns out to be a necessary 
and sufficient condition for the full correlation, where it includes information on both 
entanglement and correlation other than entanglement. Moreover, we discuss to what 
extent this criterion can give information on the entanglement of the system. The 
criterion is generic in the sense that it can be applied to mixed and pure systems. 
Also, it can be easily extended to treat the correlation of multipartite systems. 
We compare the results obtained from this criterion to those from concurrence for 
various examples and we gain agreement regarding entanglement. We believe that 
this criterion may have a wide range of potential applications in quantum information 
theory. 
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Quantum entanglement is a useful physical resource for quantum information processing. 
For instance, the predicted capabilities of a quantum computer rely crucially on entangle- 
ment and a proposed quantum cryptographic scheme achieves security by converting 
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shared entanglement into a shared secret key [2^. For both theoretical and potentially prac- 
tical reasons, it has become essential to quantify entanglement in physical systems. Thereby, 
this subject has attracted massive interest recently. Over the years, a number of different 



criteria for quantifying entanglement of physical 



systems have been developed such as nega- 



tivity of the partial transpose js], concurrence j^, 5|, correlation functions jo-si and quantum 
discord js, [lO]. However, each one of these criteria has advantages and disadvantages. 

In this Letter, we propose a simple and novel correlation criterion for the two-qubit sys- 
tem. The criterion includes information on both entanglement and correlation other than 
entanglement. Throughout this letter, for the entanglement we use the phrase " nonclassical 
correlation" , however, for the other correlation we metaphorically use " classical correlation" . 
The criterion is generic in the sense that it works for both mixed and pure systems. Also, we 
argue how one can extend it to assess the correlation of the multipartite system. Moreover, 
we discuss to what extent the information given by this criterion represents the entangle- 
ment of the system. The steps for getting the criterion can be described as follows. We 
construct the spin-correlation matrix by calculating the correlation elements of all pairs of 
the Pauli spin operators of the qubits ( i.e. ax\ (yy\ cf^z'') ■ If this matrix is non-zero, the 
proposed criterion, which can give complete information on the correlation of the bipartite, 
is the average of the non-zero elements. Nonetheless, if all the elements of the matrix are 
zero, the average will be automatically zero declaring uncorrelation. We call this criterion 
the correlation indicator and denote it by /. It is worth mentioning that, in some contexts. 



the concept of the correlation function has already been considered in the literature 
We illustrate that / 7^ is a necessary and sufficient condition for the bipartite exhibiting 
correlation. Needless to say that if entanglement exists, it will be a subset of this corre- 
lation. Furthermore, we compare the results obtained from our criterion with those from 
the concurrence C(p) for various examples. The main objective of this comparison is to 
show that when C{p) exhibits entanglement, / provides quite similar information. Besides, 
/ can include additional information on the classical correlation. It is valuable to mention 
the recent discovery that correlations other than entanglement can be responsible for the 
quantum computational efficiency of the deterministic quantum computation with one pure 



qubit [111]. Furthermore, the comparison allows us to know when the information given by 
I represents the entanglement. Precisely, we have noted that when / > 1/3 the bipartite is 
entangled, however, for /< 1/3 the bipartite may be correlated classically or nonclassically. 
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The limit 1/3 is obtained from the comparison with the entanglement of the Werner state. 
As we p roceed, we shed some light on the sudden-death of entanglement phenomenon (SDE) 
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- llTj l in terms of /. The SDE means that the entanglement in the system disappears in 
a finite interaction time, in contrast to the exponential decay of the excited states (or co- 
herence) of some systems. Despite the fact that the SDE sounds mysterious, there is no 
dynamical distinction between separable and entangled states, since quantum states can, in 
general, evolve back and forth across the boundary between separable and entangled states 



9|, |lO|. We show that / cannot detect the SDE. In this respect its behavior is quite similar 
to that of the quantum discord [ol, Q. We will describe all these findings in detail below. 

Our starting point is the generic density matrix of the bipartite, which in the standard 
basis {|ei,e2), |ei,5'2), |5'i,e2), \gi,g2)} has the form: 

/ hii bi2 6i3 bu \ 

&21 &22 &23 ^24 
6si b-i^i h 



P 



(1) 



-'31 '^32 (^33 "34 
\ ^41 ^42 ^43 ^44 

where bu + 622 + ^33 + ^44 = 1- 

Our main concern in this Letter is the evaluation of the full correlation inherited in 
the quantum system. Based on that a state of a multipartite quantum system is called 
correlated if it cannot be factorized into individual states. Mathematically, the quantum 
system described by p is correlated iff it cannot be expressed as: 



P 



Pv 



(2) 



where pj is the density matrix of the jth-party and n is the number of parties. The new 
criterion is working under the condition that Eq. ([2]) is satisfied regardless of the system 
being in the pure or in the mixed state. 

Next, we shed some light on the concurrence, which will be the benchmark for the sug- 

Eisted criterion regarding entanglement. The concurrence C{p) can be described as follows 
,&: 

C{p) =m8,^{{),y/iri- y/vl}, (3) 
where the elements Vj are the eigenvalues in decreasing order of the auxiliary matrix 



i = P{Oy ® ay)p*{ay ® CTy), 



(4) 
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where jf denotes the complex conjugation of p in the standard basis and (jy is the Pauli-spin 
matrix. We have 1 > C(p) > 0, where C{p) =0 indicates separability (zero entanglement) 
and C(p) = 1 means maximal entanglement. 

Now, we are in a position to develop the new criterion, which quantifies the correlation 
of the bipartite. For this purpose, we define the spin-correlation matrix as: 



Ixyz 



(l I I \ 

^xx ^xy ^xz 

^yx ^yy -^yz 

y I zx I zy I zz 



(5) 



where 

/.. = l(^f^^'^)|-|('^f^)(^f)l, ^,j = x,y,z. (6) 

The modulus in (|6]) is given to achieve a quantitative description for the correlation of the 
bipartite. In this regard, the elements are quite different from the covariance ones given 
in [3] . The reason for using all possible elements lij in the matrix ([5]) is that the occurrence 
of correlation (and consequently entanglement) at a time instant t necessarily implies the 
existence of at least a couple of operators a'f ■* acting on the bidimensional Hilbert spaces of 
spin 1 and 2 such that the corresponding element is non-zero. When Jjj 7^ 0, i 7^ j or la ^ 
means that the bipartite is correlated in the ij-plane or in the i-axis direction. Actually, 
there is a similarity between the structure of the matrix ([5]) and that of the moment of inertia 
tensor matrix of the 3 dimensional rigid body as well as that of the covariance matrix of a 
trivariate random vector whose probability density function is proportional to the pointwise 
density of the rigid body. For any arbitrary quantum atomic state, one can easily verify 
that crf^)\ > |(<5"i^'')(o"j^^)| for all Also, the modulus of the expectation value of 
any multiples of the Pauli spin operators always exists inside or on the surface of the Bloch 
sphere. Thereby, we have < Iij < 1, where lij = (1) means that the bipartite is 
disentangled (maximally correlated/entangled) in the ij-plane. Moving forward, we show 
that the matrix ([5]) includes all necessary information about the correlation/entanglement 
of the bipartite. The elements of the density matrix ([1]) can be categorized into three 
groups, namely, diagonal, anti-diagonal {&14, 6235 ^325 &4i}5 and the remaining elements. The 
first group contributes to the element Izz, the second one gives information on {crl^^al'^^) 
with i = x,y, while the last group provides information on {crl^^a^^^) and {o'l'^^) where 
i,j = x,y,z,i 7^ j, I = x,y and k = 1,2. In order to obtain full information on the 
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correlation of the bipartite, one has to deal with all elements in the matrix ([5]), not only some 
of them. This contrasts with the early study in ^, where the authors confined themselves 
to two elements only. We elaborate this point by giving the following example. Assume that 
we would like to assess the correlation of the of the quasi-Bell-state: 



1^) = cos6'|ei,e2) + sin6'|5fi,5f2). (7) 

If we restrict ourselves to the elements h^j-, we will find that the bipartite is disentangled, 
which is not true. Considering all these facts, the new criterion (, i.e. the correlation 
indicator /), which can give accurate information on the correlation (classical or nonclassical) 
of the bipartite, is: 

where represents the number of the non-zero elements in Ixyz- As we mentioned above, 
if lij = for all i and j, then 1 = and the bipartite is disentangled. The indicator / 
ranges between (disentangled) and 1 (maximally correlated/entangled) bipartite. It is 
obvious that / has various merits. It is a well-behaved quantity and can be applied for both 
pure and mixed systems. It is easy to be manipulated since it deals with the moments of 
the spin operators. It gives a simple means to interpret and to measure the correlation. 
Furthermore, the criterion ([8]) can be easily extended to treat the correlation/entanglement 
of the multipartite systems. 

Our next task is to check the efficiency of the criterion ([8]). This may be done by applying 
it to some examples and comparing the obtained results with those from the concurrence. 
To achieve this goal we give five examples; two of them (examples 4 and 5) are devoted to 
the SDE. 



Example (1): Assume we would like to quantify the corre^ 



three-particle Greenberger-Horne-Zeilinger state (GHZ) 18|, which has the form, e.g.,: 



ation between two qubits of the 



1"^+) = -^[\ei, 62,63) + 191,92,93)]- (9) 

The density matrix of any two qubits of the state ([9]), say 1, 2, can be evaluated by tracing 
over the third particle as: 

P = ^[|ei,e2)(e2,ei| + \9i, 92) {92, 9i\]- (10) 
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Foremost, the density matrix ( ITOj) represents a maximally mixed state. In terms of the 
concurrence, state ffTOj) is disentangled, where C(p) = 0. This is related to the definition of 
C{p), which is based on the spin flip transformation. Nevertheless, state (JTOl) is maximally 
correlated with respect to our criterion, where 1 = 1. 

Example (2): We study the entanglement of the quasi-Bell-state (JTj). For this state, the 
indicator / can be expressed as: 

/ = + + = ^[sm\2e) + 2| sin(2^)|]. (11) 

Also, its concurrence takes the form: 

C{p) = sm\2e). (12) 

It is evident from equations (fTTj) and (fT2!) that both / and C(p) have almost the same 
behavior. In this case, / gives information on the entanglement. 

Example (3). Here we shed some light on correlation of the general case of 2-Q2bit pairs. 



Also we compare the results obtained from our criterion with those given in 19|, 20 1 as well 



as with the concurrence. To this end, we use the general form of a 2-q-bit state 



19 



+ + t.r^ + ^.^.r^], (13) 



where '^C = (a^l^). We employ the terminology and the notational conventions of 19|, |20|, 



in which the individual q-bits are described by means of analogs of Pauli's spin vector 

:he second. The arrows — denote row vector, 
it has been shown that the characterization of 



operator: ~^ for the first q-bit, for 



column vector and the cross dyadic. In \19 



the 2-q-bit states produced by some source requires the experimental determination of 15 real 
parameters. Based on the knowledge of these parameters, the 2-q-bit state has been classified 
into six classes of families of locally equivalent states. Simple criteria have been stated for 
checking a given states positivity and separability through fulfilling certain inequalities. The 
technique has been employed in finding the Lewenstein-Sanpera decompositions (LS) [21 1 
for the state p. This will be our main concern in this comparison. Precisely, LS states that 
any 2-q-bit state p can be written as a mixture of a separable state psep and a non-separable 
pure state Ppure- 

P = Apsep + (A - l)ppure- (14) 
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There are many different such LS decompositions with varying values of A. Among them is 
the unique optimal decomposition, the one with the largest A value, 

~P = P& + {'^-p)P^Z /i = max{A}. (15) 

In this case /i represents the degree of separability of p; the value p = obtains only if 
p itself is a nonseparable pure state. For the separability of ( !T3ll . we restrict ourselves to 
two examples, which have been studied in detail in 19 . 
rank-2 states. 



20|, namely, Werner states and the 



For Werner states, the generic density matrix (fT3|) takes the form: 

1 ^ 

p^„ = (1 - x)prMaos + XpBell = ~ x'^ .^C^lt k-T^], (16) 

fc=l 

where the subscripts choas and Bell stand for the 2-qubit density matrices of the thermal 
and the Bell state, respectively. For the sake of clarification, in the transition from fll3p to 
flT6|) o ne has s-l- = 0, t = 0, = —x ^ e\lt k- For this state, the degree of separability 
p |l9|, concurrence C(p) and correlation indicator / for ( IT6l) can be expressed as: 

p=<(: if-3<^<3' c{p)=mBx{Q,hx-\)}, I=\x\. (17) 

1(1 -a;) ifi<x<l, 2 3 

The result of p agrees well with both the numerical findings of Lewenstein and Sanpera 



2l| and the concurrence. Moreover, it is clear that p,C(p) and / agree with each other 



for the extreme values. It should be borne in mind that I has been derived based on the 
definition Based on that, the correlation found in the interval — | < x < |in/ can 
be understood classically. Inspired by this observation we can modify /, for this case, to 
provide the nonclassical correlation in ^ < a; < 1 by subtracting 1/3, i.e. / — 1/3. However, 
the modified version needs to be normalized to unity at X — 1 5 cLS it should be. Then we get 
I = C{p). It is evident that the degree of separability p and / give identical information 
on the entanglement of pw Also, one can conclude: for any bipartite when J > 1/3 the 
correlation can be interpreted as entanglement. 
Now we give attention to states of rank 2 20|: 



Prank2 = ^[l + (cTs + XT^) ^Va{9) + ((TiTi - X(T2T2) COs{9) + Xa^T^], (18) 
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with —1 < X < 1. The generic form for states of rank 2 can be found in 19|. The degree of 
separabihty, concurrence and correlation indicator / for f|T8|) are: 



, 1 if cos(^) = 0, 

yU = < C{p) = |xcos(6')|, 

1 - |x| if cos(^) ^ 0, (19) 

/ = + |x|)| cos(^)| + |x| cos2(^)]. 

It is obvious that these three criteria provide the same information for the cases cos(^) = 
and X = cos{6) = 1. However, for x = 0,cos(6') = 1, the bipartite is separable in terms of 
fi and C{p), but / = 1/3. This confirms the above observation: the correlation exhibited 
by / is nonclassical whenever / > 1/3. For certain range of the parameters, C{p) and / 
have quite similar behavior, but this is not the case for p. This is noticeable by comparing 
different expressions in f lT9|) . As an example, when x = l,cos(0) = 1/2 we have p = 
(maximum entanglement), whereas C{p) = 1/2 and / = 5/12. 

Example (4): This example and the following one have been very often used in the lit- 



erature to establish the concept of the SDE 12l-ll7l|. Thus, we use them to discuss the 



SDK in the framework of /. We consider the evolution of the state ([7]) with the two-mode 



Jaynes-Cummings model 



ISHlTJ, which may be described by the interaction Hamiltonian: 



2 



i^, = A^(aJ\ + ai^^), (20) 

where and cr^'' {(ij) are the raising and the lowering atomic (field) operator of 

the jth-party, respectively, A is the coupling constant and we take h = 1. For the initial 



field state |0i,02), the elements of the matrix ([T]) of the two qubits take the forms 16| : 
bii = cos^ 9 cos^ T, 622 = ^33 = i cos^ 9 sin^(2T), 644 = cos^ 9 sin^ T + sin^ 9, 614 = 641 = 
i sin(2^) cos^ T, with T = tX, and the other elements are zero. The concurrence and the 
indicator / of the two qubits of this system read: 

C{p) = 2 cos'^Tcos'^ 9 max{0, Itan^l - sin^T}, (21) 
J = 1{1 - sin2(2T) cos^ 9 - [cos(2T) cos^ 9 - sin^ 9]^ 

+ 1 sin(2^) cos2 T + i cos^ 9 sm'^{2T) \ + \ sm{29) cos^ T - i cos^ 9 sin2(2T) |}. 

Comparison between (pT]) and (122|) reveals that C{p) exhibits SDE only when |tan6'| < 
sin^T, while / cannot. The reason is that / includes the classical and the nonclassical 
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correlation, while C{p) has the nonclassical one. In Figs. 1(a) and (b) we depict / and C{p) 
for two different cases 9 = /A and 7r/6, respectively. From Fig. 1(a), the behavior of the two 
criteria is almost the same, where, for this case, C(p) = cos^ T and I = |(2 + cos^ T) cos^ T. 
Fig. 1(b) shows that / has smooth oscillatory behavior without SDE, whereas C{p) non- 
smoothly becomes and stays zero for a finite interval of time. As a result of the lossless nature 
of the evolution, the original entanglement value occurs in a periodic way following each 
sudden death event. The beauty of this example is that SDE occurs without decoherence of 
the traditional type. The origin of the occurrence of SDE in C{p) is in the cut-off existing in 
the definition of the concurrence ([3]). As we proceed, one can observe that / vanishes only 
at discrete instants, i.e. T = 7r/2, declaring the disentanglement in the bipartite. One can 
verify this result by substituting (T, 0) = {it/2, 7r/6) into the density matrix of the bipartite, 
which reduces to p = \gi){gi \ |5'2)(fi'2|5 i-e. the bipartite becomes pure and separable. In 
other words, / = means that the information is locally accessible and can be obtained by 
distant independent observers without perturbing the bipartite state. The horizontal dashed 
line in Fig. 1(b) confirms our observation: for /> 1/3 the bipartite is entangled. 
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FIG. 1: Atom-atom entanglement quantified by C{p) (solid curve) and / (dashed curve) when the 
atoms are initially prepared in state ([7]) with 9 = ir/i (a) and 7r/6 (b). The vertical and horizontal 
dashed lines in (b) to show the interval of the classical correlation in / or the the segment of SDE 
in concurrence. The time scale is the vacuum Rabi period. 

Example (5): Here, we consider two two-level atoms coupled individually to two cavities 
which are initially in their vacuum states. In the framework of system-plus-environment, 
the two two-level atoms are identified as the system of interest, whereas the two cavities 
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serve as the environments. The total Hamiltonian describing this scenario is 13| : 



(23) 



+ 



where gk-, fk are couphng constants and the other notations have the same standard meaning. 
When the atoms are initially entangled with each other but not with the cavities, the solution 
of the system has been derived already in |l3| via the Born-Markov approximation. We 
focus our attention on the case that the atoms are initially in the mixed entangled state 
with 611 = a/3 > 0,622 = &33 = &23 = &32 = 1/3) ^44 = (1 ~ c^)/3 and the other elements 
are zeros. Therefore, the dynamical density matrix elements of the two atoms have the form 

m 



buit) = l^hla, b22{t) = 1(7? + ifnla), hsit) = ^(72 + 72^2^), 



(24) 



644(t) = l{l-a + Ql + Ql + QjQla), h^it) = h^it) = |7i72, 

and the other elements are zero. In using the Markov limit results and assuming the cavities 

exp(-n/2), Qi 



are similar, we have 71 = 72 = exp(— rt/2), = 1^2 = v 1 — exp(— Ft), where the T is 
the Einstein A coefficient for the two-level atoms in the cavities. The concurrence and the 
correlation indicator of this case are 



13(: 



C(p(T)) = max{0, \b23{T)\ - y/bu{T)b,^{T)}, 

I{T) = lb23iT) + 1[1 - 4622(T)] - 1 [1 - 26n(T) - 2622(T)]' 

where T = Ft. For the initial state the expressions fl2Sl) reduce to 



(25) 



C(p(0)) = -[1 - vMl^], m = ^(2a - 4a' + 7). 



(26) 



It is evident that the initial two qubits are always entangled in terms of the two criteria 
(where / > 1/3) regardless of the values of a. In Fig. 2(a) and (b) we have plotted C{p{T)) 
and /(T), respectively. Fig. 2(a) presents the well-known figure of the SDK. This figure 
shows that within the general exponential character evident in ( 12^ . disentanglement can be 
completed, i.e. C{p{T)) = 0, in a finite time while the local decoherences need an infinite 
time. For instance, when a = 1, this finite time is T = In(^^Y^). The comparison between 
figure 2(a) and 2(b) reveals when C(p(T)) exhibits entanglement, I{T) ^ declaring that 
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the bipartite is correlated. The behavior of the C{p{T)) and I{T) are quite different. For 
/(T), at T = , I{T) exhibits its maximum value as given by fj26|) . As the interaction time 
increases, /(T) goes down rapidly until local minima are achieved, which result from the 
property of the modulus " |...|" involved in the definition of the lij. After that I{T) decreases 
asymptotically to zero, i.e. I{T) cannot exhibit SDE. It would be more convenient to find 
the interaction time after which /(T) = 0. By means of (!24|) - (!25l) . one can easily prove that 
the solution of this equation is exp(— Ft) = 0. This means that /(T) needs infinite time to 
vanish. From this example and the above one, one can see that it is impossible to find sudden 




(a) (b) 

FIG. 2: Atom-atom entanglement quantified by C{p) (a) and / (b) when the atoms are initially 
prepared in the mixed entangled state with 6ii = a/3 > 0,622 = 633 = 623 = 632 = 1/3, 644 = 
(1 - a)/3. 



death of correlation/entanglement in the framework of the /. The reason for this is obvious: 
I carries information on both classical and nonclassical correlation. More illustratively, SDE 
occurs when the entanglement evolves abruptly to the classical one, however, the latter is still 
noticeable in I. It is worthwhile mentioning that this behavior, the disappearance of SDE, 
has already been addressed also via quantum discord for dissipative two-qubit dynamics 
within Markovian and non-Markovian environments jol, [lo| . Choosing initial conditions that 
manifest SDE, the authors there have compared the dynamics of entanglement with that of 
quantum discord. The analysis showed that in all cases where the entanglement, in terms 
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of the concurrence, suddenly disappears during a finite time interval, the quantum discord 
vanishes only in the asymptotic limit or at discrete instants. This reveals that quantum 
discord is more robust than the concurrence against SDE. 

It is imperative to extend the definition of the bipartite indicator / to multipartite 
systems. Indeed, in the multipartite case, apart from fully separable and fully corre- 
lated/entangled states, there also exists the notion of partial separability. Here, we focus 
our discussion on full correlation/entanglement, however, the partial one can be treated 
similarly. For multipartite systems, one has to deal with all permutations of the correlation 
elements of the spin operators. For instance, for the tripartite, one should consider the cor- 
relation elements of both two particles and three particles. However, this cannot be done so 
easily. Alternatively, one may quantify the correlation for each two qubits individually using 
dH]) and eventually carry out the average over all indicators. For fully correlated/entangled 
states, this has to be done under the condition that the correlation must be available in 
each pair, i.e. none of the individual indicators equals zero. In this approach one usually 
exploits the pairwise symmetry of multipartite system, which simplifies the calculation of 
these indicators considerably. In order to verify this technique, it is sufficient to compute 
the correlation of the three particle GHZ ffTOl) . which reads: 

/ = i[/(l,2)+j{2,3)^j(l,3)]^l^ (27) 

3 

where the indicator J^*'-'^ denotes the correlation between particle i and particle j, and we 
have used the results of the example (1). The result given by fl27|) represents nonclassical 
correlation and this agrees with the fact that the GHZ is a maximally entangled tripartite 
state. Nonetheless, this result has already been deduced in terms of the 3-tangle t^, (or 
residual tangle) 22], however, vanishes for a particular class of entangled three-qubit 
states [231]. This class is represented by the state \W): 

\W) = -^[]ei,5f2,fi'3) + \gi,e2,g3) + ]5'i,5'2,e3)]. (28) 

By applying our criterion to (125]) . we obtain / = 5/9, i.e. the W state is entangled but not 
maximally. This shows that the GHZ and the W state are locally inequivalent. 

In conclusion, we have developed a new criterion for quantifying correlation and/or en- 
tanglement of the bipartite systems. The proposed measure is simple and involves only 
calculation of average values of the spin matrices, without the need for operations on the 
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density matrix of the system nor the calculation of eigenvalues, as it happens with the 
concurrence. Also we have explained how it can be extended to treat the correlation of 
multipartite systems. The criterion is based on observable quantities, which have a clear 
physical meaning. This can give reliable physical interpretation for entanglement. Also it 
suggests a way for measuring the correlation experimentally. We have compared the results 
obtained from our criterion with those from the concurrence, and for particular examples 
with the degree of separability. There is agreement and disagreement between the behaviors 
of both criteria. Generally, there is agreement among the criteria in the entanglement do- 
main. The C (p) and I could have similar behavior when C (p) behaves smoothly. Moreover, 
from the comparison with the other criteria, we have found, when J > 1/3, the bipartite is 
entangled, however, for / < 1/3 the bipartite may be correlated classically or nonclassically. 
Also, we have shown that the indicator I cannot exhibit SDE. In this regard, the behavior 
of . . ..U. .o t.a. of d.co.d BQ. Now. .Ke <,ue.tion: . f.e.e .e,a- 

tionship between I and quantum discord? This point will be studied in a future work where 
more general Hamiltonians and initial conditions will be considered. From the analysis given 
in this Letter one can arrive at the following fact: if it is proved that a physical system is 
entangled in terms of any criterion, it will be automatically entangled in the framework of 
J, i.e. / 7^ 0, however, the reverse is not true. Finally, we believe that the criterion given in 
this Letter may have a wide range of potential applications in quantum information theory. 
For instance, the experimentalists always measure the correlation in the quantum systems 
regardless it is classical or nonclassical , since they cannot distinguish between them. In this 
sense, the new criterion could be of interest. 
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